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ABSTRACT 


We observed three massive subhalos in the Coma cluster with Suzaku. These subhalos, labeled 
“ID 1 ”, “ID 2”, and “ID 32”, were detected with a weak-lensing survey using the Subaru/Suprime- 
Cam (|Okabc ct al.l l2014al) . and are located at the projected distances of 1.4 r 50 o, 1.2 r 50 o, and 
1.6 r 5 oo from the center of the Coma cluster, respectively. The subhalo “ID 1” has a compact 
X-ray excess emission close to the center of the weak-lensing mass contour, and the gas mass to 
weak-lensing mass ratio is about 0.001. The temperature of the emission is about 3 keV, which 
is slightly lower than that of the surrounding intracluster medium (ICM) and that expected for 
the temperature vs. mass relation of clusters of galaxies. The subhalo “ID 32” shows an excess 
emission whose peak is shifted toward the opposite direction from the center of the Coma cluster. 
The gas mass to weak-lensing mass ratio is also about 0.001, which is significantly smaller than 
regular galaxy groups. The temperature of the excess is about 0.5 keV and significantly lower 
than that of the surrounding ICM and far from the temperature vs. mass relation of clusters. 
However, there is no significant excess X-ray emission in the “ID 2” subhalo. Assuming an infall 
velocity of about 2000 km s _1 , at the border of the excess X-ray emission, the ram pressures 
for “ID 1” and “ID 32” are comparable to the gravitational restoring force per area. We also 
studied the effect of the Kelvin-Helmholtz instability to strip the gas. Although we found X-ray 
clumps associated with the weak-lensing subhalos, their X-ray luminosities are much lower than 
the total ICM luminosity in the cluster outskirts. 


Subject headings: galaxies:clusters:individual(Coma Cluster) -X-rays:galaxies:clusters 


1. INTRODUCTION 

Galaxy clusters are the largest self-gravitating 
bound systems in the Universe and composed of 
thousands of galaxies, the intracluster medium 
(ICM), and the dark matter. The ICM cov¬ 
ers the total mass range of roughly 10 13-14 M 0 
and is bound to the potential of the dark mat¬ 
ter halo which covers that of 10 14-15 Mq. There- 
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fore, the gravity of the dark matter halo plays the 
most important role in cluster evolution and struc¬ 
ture formation. Numerical simulations with the 
Cold Dark Matter (CDM) model predict that the 
galaxy clusters form through merger or accretion 
of the smaller system like galaxy groups. Since the 
dynamical time scale of galaxy clusters is compa¬ 
rable to the Hubble time, outskirts of the galaxy 
clusters still maintain the evolution effects via the 
accretion of the substructures. The central region 
of these accreting objects is expected to be sur¬ 
vived until recent days as subhalos in the cluster 
host halo. 


The mass distribution of subhalos provides 
with information on the mass as sembly of 

( 2014a ) surveyed 


us 

galaxy 

and measured the mass of subhalos in the Coma 


cluster. Okabe et al 
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cluster using weak-lensing observations with the 
Subaru/Suprime-Cam. Thanks to the large appar¬ 
ent size, they detected 32 cluster subhalos whose 
mass range is ~ 2 — 50 x 10 12 h~ 1 M Q . They 
first confirmed that the subhalo mass function, 
d n/d In M su b, is well represented with a single 
powerlaw or a Schechter function. The best-fit in¬ 
dexes of each model are ~ 1, which agree well with 
CDM model prediction on sub-scale of cluster. 
Stacked signals of subhalos were well represented 
with a sharply truncated Navarro-Frenk-White 


(NFW) mass model (INavarro et al.l 119951 ) as ex¬ 


pected from a tidal destruction model. For three 
most massive subhalos whose mass are higher than 
~ 1 x 10 13 they measured mass and trun¬ 

cation radius of each subhalo. One of them is a 
famous substructure of the Coma cluster around 
the NGC 4839 group. 

If subhalos in the cluster outskirts still pos¬ 
sess some amount of their hot gas and are not 
excluded in X-ray analysis, the derived ICM den¬ 
sity would be overestimated. The recent Suzaku 
observations reported that the entropy of the 
ICM, which is a useful parameter for the ther¬ 
modynamical history, shows flatter profiles be- 


yond rr,nn than expectations from pure gravita- 
tional heating (George et al. 2009t Rciprich et al. 

2009: Bautz et a 

. 20091: Hoshino et al. 

2010 


Kawaharada et al 

201(ir Simionescu et al. 

2011 

Urban et al.l201ll 

Akamatsu et al. 2011; Wa 

ker et al. 

2012al: Akamatsu et al.ll2012l 

Walker et al. 20124 

Sato et al.l 2012- Walker et al 

2013; Ichikawa et al. 

2013: Simionescu et al. 2013 

Illrban et al.l 2011 

; 

Sato et al. 2014); Okabe et al. 20144. Since the 


gas fraction in the Perseus cluster assumed hy¬ 
drostatic e 2 uihbnum_exceeds_the cosmic baryon 


fraction, Simionescu et al. ( 201 lh claimed the ef¬ 


fect of the gas dumpiness, although there is no 
significant excess X-ray sou rces in the outskir ts 


with a Chandra observation (Urban ct al.ll2014 ). 


Weak-lensing mass estimation of cluster main 
halos are complementary to X-ray observations. 
In the cluster outskirts, the hydrostatic mass with 
Suzaku is significantly lower tha n the weak -l ensing 


mass with the Subaru teles c ope |K awaharada_et_aL 
2010: Ichikawa et al.l 12013 : Mochizuki et al. 201.4 


Okabe et al.l 2014bT . Okabe et al. ( 2014b ) dis¬ 

cussed that the bivariate scaling functions of the 
electron density and temperature indicate that 
entropy flattening of the outskirts of the galaxy 


clusters caused by the steepening of temperature 
profiles. Deviations from h ydrostatic e quilib rium 
have been disc u ssed in iKawaharada et al.1 ( 2010 ): 
Ichikaw^etaL (120131) : iMochizuki et ahl (2015 ); 
Okabe et al.l ( 2014bl) . Hoshino et al.l ( 2010h pro- 
posed an another idea that electron temperature 
is lower than the ion temperature in these regions, 
since heating the electrons takes a longer time 
than that of the ion after accretion shocks and 
mergers. 

With weak-lensing mass measurements of sub¬ 
halos, we can search X-ray clumps associated with 
these subhalos efficiently with X-ray observations. 
In this paper, we describe the X-ray properties of 
three massive subhalos, whose mass is greater than 
~ 9 x 10 12 h~ 1 M Q , detected Subaru weak-lensing 


observations of the Coma cluster (Okabe et al 


2014a ). Excluding the NGC 4839 subgroup, we 
observed two of the most massive subhalos with 
Suzaku. We also observed a smaller subhalo, 
whose total mass is ~ 9 x 10 12 h~ 1 M@. We sum¬ 
marize the observations and date preparation in 
section [2] Section l3Tl shows the X-ray images and 
surface profiles of each subhalo. In section EPl we 
summarize the spectral fitting and that results. 

We compare of the X-ray properties with those of 
other galaxy groups, and discuss the effect of ram 
pressure and gas dumpiness in section |4] 

We use = 0.27, 12 a = 0.73, and H 0 =70 km s -1 Mpc -1 
in this paper. At the redshift of the cluster, 2 = 

0.0231 ( Struble fe Rood 1999h . I aremin corre¬ 
sponds to 28.9 kpc . The solar abundance table 
is given bv iLodders ( 2003 ). The errors are in the 
68% confidence region for the single parameter of 
interest. 


2 . 


OBSERVATION AND DATA REDUC¬ 
TION 


In the Coma cluster, Okabe et al. (2014a}) de¬ 
tected three massive subhalos whose mass are 
higher than 10 13 M 0 ft -1 , and which are labeled 
as “ID 1”, “ID 9”, and “ID 32”. The weak- 
lensing signal of these three subhalos were well 
represented by a truncated NFW model. Since 
the “ID 9” is associated with a halo of the fa¬ 
mous subgroup around NGC 4839, which has 
be en already observ e d wit h Suzaku and reported 
bv lAkamatsu et al ( 20131) . we observed “ID 1”, 
“ID 32”, and a southern offset region of “ID 32” 
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Table 1 

Properties of the Coma cluster subhalos. 


ID a 


M c 

r? 

(R.A., decl.) e 

ATf 

10^ cm - 2 

Distance 9 


10 12 h- 1 Mq 

10 12 h.- 1 M 0 

arcmin 

J2000.0 

Arcmin/r 50 o 

1 

15.42 ± 2.79 

14 26' 1 ’ 2 ' 37 

1 ^' ZD -2.53 —5.55 

3.86t°;“ 

12 h 55 m 34®5, +27°3l'33.7" 

8.6 

61.8/1.42 

2 

8.79 ± 4.69 

- 

- 

12 h 56 m 03 s 8, +27°47'20.8" 

8.7 

51.6/1.18 

32 

45.95 ± 7.57 

47 75+ 5 - 81 

• 1 -5.81-13.42 

9 21" 1 " 0 ' 74 

13 h 01 m 41?0, +29°03'14.4" 

9.5 

71.2/1.63 


a The name of subhalos dOkabe et al.|[2014aD . 

b The projected weak lens mass of the subhalos dOkabe et al.|[2014af) . 
c The best-fit mass with truncated NFW model dOkabe et al.|[2014a|) . 
d The truncation radius derived from dOkabe et al.|[2014aD . 

e For “ID 1” and “ID 2”, the center of the subhalos determined from the mass contour. The center of “ID32”, however, 
is derived from the weak-lensing signal peak since the difference between weak-lensing signal peak and mass contour is 
significantly smaller than the Suzaku point spread function. 

f The Galactic hydrogen column density (|Kalberla et aOl2005li . 

g The distance from the X-ray peak of the Coma cluster center (12 h 59 m 44®81, 27°56 / 4{P92). 


Table 2 

Suzaku OBSERVATION LOGS FOR COMA CLUSTER SUBHALOS. 


Field name 

Sequence 

Number 

Date-Obs. a 

(R.A., decl.) b 

J2000.0 

Exposure 0 

ksec 

ID 1 

808022010 

2013-06-10T14:12:00 

12 h 55” 

“28?0, 27°3l'00"l 

18.4 

ID 2 

808021010 

2013-06-10TOO:39:15 

12 h 55” 

“55? 3, 27° 45' 17" 6 

23.7 

ID 32 

808018010 

2013-06-08T09:04:42 

13 h 01” 

a 36!l, 29°0l'40"8 

26.4 

ID 32 BGD 

808019010 

2013-06-09T01:38:51 

13 h 01” 

a 00!7, 28°45'46"4 

20.2 


a Start date of observation, written in the DATE-OBS keyword of the event FITS files. 

b Average pointing direction of the XIS, written in the RA.NOM and DEC_NOM keywords of 
the event FITS files. 

c Exposure time after screening. 
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as a background (hereafter “ID 32 BGD”) with 
Suzaku. We also observed “ID 2” subhalo ,which 
is associated with the NGC 4816 group and the 
total mass reach to 9xl0 12 Mq fa -1 , with Suzaku. 
The mass, truncation radius (hereafter r t ), and co¬ 
ordinates of each subhalo are summarized in table 
[Hand observational logs with Suzaku are shown in 
table [2] The three subhalos, ”ID 1”, ”ID 2”, and 
”ID 32”, are located at the projected distances 
of 1.4 f 5 oo, 1.2 7 - 500 , and 1.6 r $oo from the X-ray 
peak of the Coma cluster, respectively. In figure 
[TJ we overlaid the field of views (FOVs) of Suzaku 
pointings and the contours of mass m ap derived 
from weak-lensing ( Okabe et al.ll2014ah on the X- 
ray image with ROSAT. We also used four Suzaku 
pointings beyond 2.5 rsoo of the Coma cluster to 
study the X-ray background emission. The details 
are described in Appendix [A] 

In this study, we used only XIS data. The XIS 
instrument consists of three sets of X-ray CCDs 
(XIS 0, 1, and 3). XIS 1 is a back-illuminated (BI) 
sensor, while XIS 0 and 3 are front-illuminated 
(FI). The instruments were operated the normal 
clocking mode (8 s exposure per frame). The data 
were reprocessed the standard screening criteria 
0 using HEAsoft 6.15. We also performed event 
screening with the cosmic-ray cut-off rigidity COR 
> 6 GV, and the Earth rim ELEVATION >10°. 
We generated Ancillary Resp onse F i les ( A RFs) b y 


“xissimarfgen” Ftools task ( Ishisaki ct al. 20071), 


assumed a uniform sky of 20' radius. The ef¬ 
fect of degrading energy resolution by radiation 
damage was included in the redistribution matrix 
files by “xisrmfgen” Ftools task. We employed 
the night Earth database generated by the “xis- 
nxbgen” Ftools task for the same detector area to 
subtract the non-X-ray background (NXB). 


3. ANALYSIS AND RESULTS 


3.1. X-ray images and Surface brightness 
profiles 


vignetting effect using a flat image at 1 ke\0. 
We created surface brightness profiles from these 
images of individual subhalos along the direction 
to the center of the Coma cluster, (R.A., Deck) 
= (12 h 59 m 44!81,27°56'49"92). The resultant sur¬ 
face brightness profiles are shown in figure [3] 

As shown in figure [2] (a) and figure [3] (a), an 
excess emission is seen around the center of the 
”ID 1” mass contour. The surface brightness pro¬ 
file shows that most of the excess emission is con¬ 
fined within ~ 2'(0.6 r t ) from the mass center. In 
contrast, the ”ID 2” subhalo does not show any 
excess emission: the brightness profile gradually 
increases toward the Coma cluster center as shown 
in figure [3] (b). 

In figurc[3](c), the ”ID 32” subhalo shows an ex¬ 
cess emission whose peak is not located on the cen¬ 
ter of the mass contour, but shifts about 3' (0.3 r t ) 
away from the subhalo center toward the northern 
part or the opposite direction from the Coma clus¬ 
ter center. The excess emission is extended to at 
least 5'(0.6 r t ) from the mass center. Because of 
the asymmetrical profile for the “ID 32”, we also 
made projections of the surface brightness of the 
northern and southern parts from the mass cen¬ 
ter of the subhalo along the perpendicular direc¬ 
tion against in figure [3] (c). The resultant projec¬ 
tions of the surface brightness are shown in figure[3] 
(d). The excess emission at the northern direction 
extends out to 5-6' (0.6 r t ) toward the east and 
west directions. The peak of the X-ray emission 
is located in the northern part and shifts about 
3' (0.3 r t ) away from the subhalo center toward 
west direction. In contrast, the brightness of most 
of the southern part is consistent with that of the 
background region, although there are two peaks 
at ~ 6' offsets of east and west directions. We 
extracted spectra around these two excess. Since 
their spectra are relatively hard and fitted with a 
power-law model, the southern peaks are possibly 
caused by background point sources. 


In figure [2] we present combined XIS images 
of subhalos in an 0.5-2.0 keV energy band. Here, 
the difference in the exposure times are corrected 
with exposure map generated by “xisexpmapgen” 
Ftools task. In addition, we also corrected the 


http: //heasarc.nasa.gov/docs/Suzaku/processing/criteria_xis.html 


3.2. Spectral fitting 

We extracted spectra over subhalo regions and 
background regions as shown in figure 0] The 
regions around point sources brighter than 1 x 
10 13 erg s -1 cm -2 in 2.0-10.0 keV were excluded 
from the spectral analysis. Since the mass con- 


2 http: //heasarc.gsfc.nasa.gov/docs/suzaku/analysis/expomap.html 
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Fig. 1. The weak-lensing mass contour map in linear scale from Subaru overlaid on the X-ray image taken 
by ROSAT All Sky Survey in 0.1-2.4 keV band. Here, the exposure time was corrected and instrumental 
backgrounds were subtracted. The X-ray image was smoothed by a Gaussian of a ~ 2'. The numbers in 
color bar is in units of counts s” 1 pixel” 1 . The FOVs of the Suzaku pointings are plotted with boxes (yellow). 
The solid (white) circles ind icate r^nn (~ 44') a nd r v i r ~ rgs (~ 97'), respectively, which were derived with 
weak-lening observa tions by Okabe_etal. ( 2010h . The dashed (magenta) box shows the field of the Subaru 
weak-lcning survey (lOkabe et al.1 2014a ) and the mass distribution derived from this survey was overlaid in 
contours. The crosses (green) show the positions of representative galaxies, NGC 4807, the NGC 4816 group, 
and IC 4088, for the subhalos “ID 1”, “ID 2” and “ID 32”, respectively. (A color version of this figure is 
available in the online journal.) 


5 

















(a) 


(b) 


27.800 


= 27.600 


27.400 


29.000 


28.600 


194.200 


194.000 


193.800 


193.600 




Right ascension 

IS 1 2 


Right ascension 

IS ii.b 


Fig. 2.— The weak-lensing mass contour map in linear scale overlaid on the NXB subtracted XIS images in 
0.5-2.0 keV. (a) around “ID 1” and “ID 2”, and (b) around “ID 32” in the 0.5-2.0 keV energy band. Here, 
the difference in the exposure times and the vignetting effect at 1 keV were corrected. The images were 
smoothed by a Gaussian of a =24 pixels « 25". The numbers below the color bars have units o f counts Ms -1 
pixel -1 . The contours (white) show the mass map derived from weak-lcning bv lOkabe et al. ( 2014 a). The 
direction of the Coma cluster center from the center of each subhalo is shown by the arrow. The crosses are 
the center of each subhalo as summarized in table [Q For ”ID 1” and ”ID 32”, the solid (yellow) circles show 
the truncation radii. The diamond (cyan) mark in the ”ID 32” subhalo corresponds to the X-ray peak of 
the excess emission. When creating the surface brightness profiles, the regions in the dashed (white) circles 
and boxes were excluded. (A color version of this figure is available in the online journal.) 
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Fig. 3.— The projected surface brightness profiles of (a) ”ID 1” subhalo, (b) ”ID 2” subhalo, and (c) ”ID 32” 
subhalo along the direction to the outskirts from the Coma cluster center extracted from the combined XIS 
FOV images in the 0.5-2.0 keV energy band. The side of negative number corresponds to the direction to 
the center of the Coma cluster. Here, the NXB were subtracted and the effect of the the exposure time 
and vignetting were corrected. The minus X-axis corresponds to the direction of the Coma cluster center. 
The dashed vatical lines indicate the center of each subhalo and the dotted vatical lines show the truncation 
radius of the subhalo “ID 1” and “ID 32”. (d)The projection direction is perpendicular against (c) . The 
sides of negative and positive number correspond to the direction to east and west. 
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tours of the subhalo “ID 1” in linear scale are 
elongated, we extracted spectra over an elliptical 
region, whose semiminor and semimajor axes are 
1.0 and 1.6 rt, respectively. Here, we excluded a 
circular region aroun d a background galaxy group 
( Okabe et al.ll2014al) . plotted as a dashed circle in 
figure [4] To study background emissions includ¬ 
ing the ICM contribution, we extracted spectra 
over an square region (hereafter ”ID 1 BGD”), 
excluding the elliptical region. For the “ID 32” 
subhalo, we extracted spectra over two semicir¬ 
cular regions (hereafter ’south 1 or ’north 1 regions) 
of the subhalo out to the truncation radius. The 
background spectra for the “ID 32” were extracted 
from the FOV of the ”ID 32 BGD” observation. 


The spectral fitting was carried out using 
XSPEC 12.8.lg and the extended C-statistic esti¬ 
mator. The spectra were binned to have at least 
one count per channel. We used the energy ranges 
of 0.5-7.0 keV and 0.7-7.0 keV for the BI and 
FI detectors, respectively. We excluded the en¬ 
ergy band around the Si-K edge (1.82-1.84 keV) 
because its response was not modeled correctly. 


We assumed that the X-ray emissions from 
the ”ID 1 BGD” and ”ID 32 BGD” consist of 
the Galactic emissions from the Local Hot bub¬ 
ble (LHB) and the Milky Way Halo (MWH), 
the Cosmic X-ray background (CXB), and the 
ICM of the Coma cluster. The LHB and MWH 
were m odeled with a ther mal plasma model ( apec 
model; ISmith et all 120011 ) without and with the 
Galactic absorption, apecLHB and phabscAL x 
apecMWH, respectively. Here, the phabsc al indi¬ 
cates the Galactic absorption using phabs model 
in the XSPEC package, and the column den¬ 
sity of each subhalo direction were summarized 
in table [T] The temperature of the apecLHB was 
fixed at 0.1 keV, while its normalization was a 
free parameter. We also allowed the tempera¬ 
ture and normalization of the apecMWH to vary. 
For the LHB and MWH, the redshift and abun¬ 
dance were fixed to be 0 and 1 solar, respec¬ 
tively. The CXB emission was described an ab¬ 
sorbed power-law model with a photon of index 
r = 1.4, phabs gal* power—lawcxB- The normal¬ 
ization of power — lawcxB were estimated beyond 
the virial radius of the Coma cluster as described 
in appendix [Al 


The ICM emission was modeled with a ther¬ 
mal plasma model with the Galactic absorption, 


phabsQAB x apecicM ■ The temperature, normal¬ 
ization, and abundance were allowed to vary, ex¬ 
cept for the abundance for ”ID 32 BGD” which 
was fixed at 0.2 solar. The redshift of the ICM 
component was fixed at the value of the Coma 
cluster, z = 0.0231. Thus, we used the fol¬ 
lowing model formula for the spectra for the 
”ID 1 BGD” and ”ID 32 BGD”; constant x 
(apecLHB T phabsQAB x (apec M wh + apecicM + power - Iowcxb)) ■ 
Here, the constant is a normalization parameter 
for the difference in relative normalizations among 
XIS detectors. 


For the subhalo regions, we added an extra 
thermal plasma model, phabsc al x apec su bhaio- 
We allowed the temperature and normalization 
of apec su bhaio of the ”ID 1” and ”ID 32” to 
vary. The abundance for the ”ID 1“ was a 
free parameter but that for the ”ID 32” were 
fixed at 0.2 solar, since we can not constrain 
the abundance. Even if the abundance was 
fixed to be 0.1 or 0.3 solar, the results did not 
change. The redshift of the apec su bhaio was also 
fixed at the value of the Coma cluster. Thus, 
we modeled the spectra extracted from sub¬ 
halo regions as following formula; constant x 
(apecLHB + phabs gal x (apec M wh + apecicM + apec su b haio + power 
We finally fitted the spectra extracted for each 
subhalo region and corresponding background 
region simultaneously, assuming that the X-ray 
background components have the same surface 
brightness, temperature and abundances. Here, 
relative normalizations of three XIS detectors were 
allowed to vary. 


The fitting results are summarized in table [3] 
Figure [5] shows the best fit spectra of the back¬ 
ground regions. The ICM temperatures for the 
“ID 1” and “ID 32” regions, 4.33±g;|? keV and 
5.19iJ;gg keV, respectively, are consistent with 
previo us results of sout h west and northwest direc¬ 
tions ( Simioncscu et al • 5HJ at similar distance 
from the cluster center, respectively. Although the 
error range is fairly large, the ICM abundance 
of the “ID 1” was als o consistent with that in 
Simionescu et al.l ( 20131) . 


The temperatures and normalizations of the 
Galactic components are consistent between the 
“ID 1” and “ID 32”. The normalization of the 
LHB derived from the ”ID 1” and ”ID 32” spec¬ 
tral fits is consistent with that for the region at 
110' — 130' (Appendix). However, the tempera- 
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Fig. 4.— The same images as figure [2] The elliptical region (yellow) in the left panel and semicircular 
regions (yellow) in the right panel are used in spectral extraction for the ”ID 1” subhalo, ’south 1 and ’north 1 
regions of the ”ID 32” subhalo, respectively. The solid box (yellow) in the left panel excluding the elliptical 
region and that in the right panel are used in the background analysis (”ID 1 BGD” and ”ID 32 BGD”, 
respectively). The regions within dashed box and circles (white) are excluded from these spectral analysis. 
(A color version of this figure is available in the online journal.) 


ture and normalization of MWH are significantly 
different from those derived for the 110' — 130' re¬ 
gion. When we use the temperature and normal¬ 
izations of the the Galactic components obtained 
from 110' — 130' region, and fitted the spectra of 
the subhalo regions, the temperature and normal¬ 
ization of the ”ID 1” and the ’north’ region of the 
”ID 32” did not change within statistical errors, 
although the temperature of the ’south’ region of 
”ID 32” decreased to about 0.1 keV. Considering 
the possible spatial variation of the Galactic com¬ 
ponents, we adopted the results of the simultane¬ 
ous fits using ”ID 1 BGD” and ”ID 32 BGD”. 

3.3. Fitting results of the subhalo compo¬ 
nents 

Figure [5] also shows the best fit spectra of the 
subhalo regions. The fitting results of the subhalo 
components are summarized in tabled The spec¬ 
tra for the subhalo and background regions were 
well-represented with our model formula. Without 
the subhalo components, the data-to-model ratios 
show significant excess. 


The temperatures of each subhalo component 
are lower than that of the surrounding ICM. 
For the “ID 1”, the temperature is 2.711 q 59 
keV and is cooler than the ICM component, 
4 .33^g'37 keV. The temperatures of ‘north’ and 
‘south’ of the “ID 32” are about 0.55ito i 3 keV 
and 0.29±g o| keV, respectively. These values are 
significantly lower than the surrounding ICM tem¬ 
perature, 5.19lJ’gg keV. The abundance of the 
subhalo “ID 1” component is consistent with that 
of the surrounding ICM. The luminosity of the 
“ID 1” component is about 2 x 10 41 erg s _1 at 0.5- 
2.0 keV energy range and those of the ’north 1 and 
’south 1 regions of “ID 32” are ~2x 10 41 erg s _1 
and ~ 7 x 10 40 erg s _1 at the same energy range, 
respectively. 


3.4. The representative background struc¬ 
ture fitting 


Background galaxy groups are located within 
the “I D 1” and “ID 3 2” subhalo regio ns at z = 
0.418 ( Wen et ahll2009l) and z = 0.189 (Hao et al 


2010h . respectively ( Okabe et al. 2014a ). If the 
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Fig. 5.— The spectra of (a) “ID 1 BGD”, (b) the subhalo region spectra of “ID 1”, (c) “ID 32 BGD”, (d) 
and (e) ’North’ and ’south’ regions of the subhalo “ID 32”, respectively. The spectra were rebinned here for 
display purpose only. Upper panels show the NXB subtracted XIS 1 spectra (black crosses). The subhalo 
component is plotted by (red) bold line. The ICM, CXB, LHB, and MWH components are indicated by 
(blue) dotted, (magenta) dash-dotted, (cyan) thin, and (green) dashed lines, respectively, and the gray solid 
line indicates the sum of these background components. The lower panels show the data-to-model ratios. 
The gray diamonds in the lower panels in (b), (d) and (e) show those without the subhalo component. 
In the lower panels in (b), (d) and (e), the gray diamonds show that data-to-model ratio without subhalo 
component, the ratio indicates that the subhalo component can not be negligible. (A color version of this 
figure is available in the online journal.) 
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Table 3 

The fitting spectral results. 


Background components 

Field i 

name 

Normans 

kTMWH 

^ orrn MW H 

kTicM 

ZlCM 

N orm j c M 




keV 


keV 

Solar 


ID 1 BGD 

18- lit 1.9 

0 60' 1 ’ 0 ' 11 
u.ou _ 0 09 

0 41 + 0- 08 

0.07 

4 qq+0 - 60 

0 n +°- 17 

U '- LJ -- 0.11 

Q Q'7+0.57 

_0.65 

ID 32 BGD 

! 8 .l±i;= 

0 fil "*" 0 ' 06 
u,Oi _0.05 

0 to - * -0 ' 05 

U.OU_Q 05 

5 19+ 1 - 04 

- Li7 —0.83 

0 . 2 (fixed) 

9 a.q+ 0.20 

4.00 o i 9 

Subhalo components 



redshift 

kT 

Z 

Norm a 

Luminosity^ 

Flux 5 




keV 

solar 


10 41 erg s 1 

10” 13 erg s 1 cm ” 2 

“ID 1” 


0.0231 

9 ij-i -f-u.yy 
‘ 1 — 0.59 

f) 1 0+0.30 

u.io_Q 13 

c 79 - 1 - 1.32 
' z -1.32 

-1 '7Q"l -U-dy 
( O-0.33 

1 48" t ’ u- ^ 4 

l-^o_0.29 

“ID 32’ 

’ north 

0.0231 

0.55t ° n ° 7 

0 . 2 (fixed) 

2-80± 4 ;™ 

1 76 + 017 
° —0.40 

1 44+0-14 
-0.32 

“ID 32’ 

’ south 

0.0231 

o-29tS:i? 

0 . 2 (fixed) 

2.06ti;“ 

0 67 +0 ' 29 

uu ‘ -0.48 

0 54 +0-21 

“ID 1” 


0.418 

O 09+1'^ 

0.00_o.82 

0 09" t u 

u.uu _ 0 09 

9.834;.?^ 

(8.38+;;^) X 10 2 

1 C ) 0 +^' 2A 

l.OU — o 15 

“ID 32’ 

5 north 

0.189 

0-90to “3 

0 . 2 (fixed) 

0 n +0.51 
o.Ol-Q 51 

(1.74±JJ; 23 ) x 10 2 

1 66"*" 0-35 

1 .OO -0 20 

“ID 32’ 

’ south 

0.189 

0-681S.41 

0 . 2 (fixed) 

1 no - * -3-65 

j..u»_o 41 

(5.24±|;||) x 10 1 

0 49 + 0 - 28 

W- ^ C7 —0.32 


a The normalization of the apec components divided the solid angle, lA. assuming a uniform sky of 20 / radius, 
Norm = f n e nndV / [ 4 - 7 T (1 + z) 2 D£] / Q u Xl0 -17 cm -5 4007T arcmin” 2 . 

b The energy range is 0.5—2.0 keV. 


weak-lensing signals were mostly caused by the 
corresponding background galaxy groups, the 
virial mass, M V - IT would be 1 -3xl0 15 M 0 for 
’ ’ID 1” and several times 10 14 M 0 for ”ID 32” 
( Okabe et ah! 2014a). For such massive clusters, 
we expect that the ICM temperatures and ICM 
luminosities exceed 10 keV and 10 45 erg s _1 for 
”ID 1”, respectively, and several keV and sev¬ 
eral times 10 44 erg s _1 for ”ID 32”, respectively. 
Therefore, we refitted the Suzaku spectra using 
redshifts of the background galaxy groups. As 
a result, the temperature increased to ~ 4 keV 
and 0.8 keV for ”ID 1” and the ’north’ region of 
”ID 32”, respectively, and the X-ray luminosity 
became ~ 10 44 erg s _1 and 2xl0 43 erg s -1 in the 
0.5-2.0 keV range, respectively. These temper¬ 
atures and X-ray luminosities are far below the 
expected values assuming that the weak-lensing 
signals come from the background galaxy groups. 
Thus, it is unlikely that the excess emissions and 
weak-lensing signals come from the background 
galaxy groups, and fairly are associated with the 
Coma cluster. 


3.5. The gas mass estimation 

To estimate the gas mass of each subhalo, we 
first approximated spherical symmetry and as¬ 
sumed constant density up to 2.5 arcmin and 6 
arcmin for “ID 1” and “ID 32”, respectively , 


which correspond to ~ 0.6 rt for each subhalos, 
since beyond these radii, the excess emission was 
not detected. The best-fit normalization in the 
spectral fitting leads us to derived average elec¬ 
tron density within the extracted region of each 
subhalo. The resultant average electron densities 
of “ID 1”, ’south 1 and ’north 1 regions of “ID 32” 
are (5.3 ± 0.6) x 10" 4 cm" 3 , (2.0 ± 0.8) x 1CT 4 
cm -3 , and (2.4 ± 0.4) x 10 -4 cm -3 , respectively. 
Integrating the electron densities out to 0.6 rt, 
the derived gas mass are (2.1 ± 0.2) x 10 10 M 0 , 
(6.8 ± 1.2) x 10 10 Mq, and (5.8 ± 2.2) x 10 10 M 0 
for ”ID 1” and the ’north 1 and ’south 1 regions of 
“ID 32”, respectively. 

We also estimated the electron density pro¬ 
files by deprojecting radial profiles of the surface 
brightness centered on the center of the mass con¬ 
tour and the X-ray peak of ”ID 1” and ”ID 32” 
subhalos, respectively. We integrated the electron 
density profiles out to 0.6 rt, within which the X- 
ray emissions of each subhalo are detected. Then, 
the derived gas mass of ”ID 1” and ”ID 32” are 
(1.5 ± 0.1) x 1O 1O M 0 and (9.5 ± 0.5) x 1O 1O M 0 , 
respectively. Comparing these values with those 
assuming constant electron density, the system¬ 
atic uncertainties in the gas mass would be about 
a factor of 2-3. 

Using the radial profiles of the electron den¬ 
sity profile and the temperatures of ”ID 1” and 
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the ’north 1 region of ”ID 32” derived from the 
spectral fits, radial profiles of the thermal gas 
pressure, P = n e kT out to ~ 0.6 r*, were esti¬ 
mated and plotted in figure [ 6 ] Since the surface 
brightness and temperature of the ICM compo¬ 
nent surrounding ”ID 1” and ”ID 32” are close 
to those of the southwest direction and the az¬ 
imuthal averag es excluding the s outhw est direc¬ 
tion derived bv ISimioncsc u et al.l ([20131 1. we com- 
pared our pressure prof iles with those derived by 
Simioncscu ct al.] ( 20131) . The thermal pressure at 
0.5 rt of ”ID 1” is slightly higher than that of the 
southwest direction. That in 0.3-0.5 rt of ”ID 32” 
is slightly below the ICM pressure. 

We also estimated the gas mass assuming that 
excess emission come from background galaxy 
groups located within the “ID 1” and “ID 32” 
subhalo regions. The calculation methods are the 
same as mentioned in the previous paragraphs. 
Assuming a flat density profiles, the total gas mass 
for “ID 1” and “ID 32” are (1.8 ±0.2) x 10 13 M 0 
and (1.3±0.8) x 10 13 M 0 , respectively. By depro- 
jecting the surface profiles, we estimated the gas 
mass by integrating the calculated density profiles. 
The resultant gas masses are (1.3±0.1) x 10 13 M 0 
and (1.3T0.1) x 10 13 M 0 for “ID 1” and “ID 32”, 
respectively. 

4. DISCUSSION 

We observed three massive subhalos detected 
by the weak-lening survey with Subaru, “ID 1”, 
“ID 2” and “ID 32”, which are located on the pro¬ 
jected distances of 1.4 rsoo, 1-2 r 500 , and 1.6 r 500 
from the center of the Coma cluster, respectively, 
with Suzaku. The excess X-ray emission have been 
detected from “ID 1” and “ID 32”, while “ID 2” 
has no significant excess emission. Temperatures 
of these subhalos were lower than that of the sur¬ 
rounding ICM. In section 14.11 we compare the 
Mg as -Mtotai and kT-M to t a i relations between sub¬ 
halos and regular galaxy groups. We estimate the 
ram pressure and mass loss rate caused by the 
Kelvin-Helm holtz instability i n sect ion 14.21 As 
discussed in ISimionescu et ah ( 2011 b the clump¬ 
ing in the galaxy clusters lead us to overestimate 
the electron density. In section 14.31 we study the 
effect of the subahlos on the density and temper¬ 
ature measurements. 


4.1. Comparison of the gas mass fraction 
and temperature with other clusters 

The gas mass and weak-lensing mass of ”ID 1” 
and ”ID 32” are plotted in figure [7] (a). The gas 
mass fraction, or the gas mass to weak-lensing 
mass ratio, of these two subhalos are about 0 . 001 . 
Here, we used the gas mass derived from integrat¬ 
ing the radial profiles of electron density. We com¬ 
pared these gas mass fraction of the subhalos with 
the gas mass to hydrostatic mass ratios of clus¬ 
ters and groups of galaxies. Since the mean den¬ 
sity within the truncation radius of subhalo “ID 1” 
and “ID 32” are about 6600 times and 27000 times 
higher than the critical density of the Universe, re¬ 
spectively, we calculated gas mass and hydrostatic 
mass at a radius within which each average density 
is the same as each subhalo’s overdensity , using 
Chandra results bv IVikhlinin et ah ( 2006 1. The 
gas mass fractions of these clusters correlate well 
with the hydrostatic mass, and are about 0 . 02 - 0 . 1 , 
which are about 1-2 orders of magnitudes higher 
than those of subhalos. If these subhalos had been 
regular groups before infalling onto the Coma clus¬ 
ter, they should have lost most of their gas. 

In figured] (b), we also compared the relation 
of the temperature and total mass ( kT-M tot& \ re¬ 
lation) of the subhalos and clusters. Here, we also 
calculated the temperatur e and hydrostatic mass 


of the Chandra clusters in Vikhlinin et al. ( 2006 1 


at the radius of the same over densities, or r66oo f° r 
”ID 1” and 727000 f° r ”ID 32”. The temperature 
of ”ID 1” is slightly lower than expected by the 
kT — Mtotai relation of clusters. For “ID 32”, the 
temperature of the subhalo is several times lower 
than those of the kT — M to tai relation of clusters. 

4.2. The effect of ram pressure stripping 

The observed very low gas mass fraction and 
morphologies of the excess emission indicate that 
the gas in the infalling subhalos has been stripped 
via ram pressure of the surrounding ICM. The sub¬ 
halo would be unable to hold the interstellar ma¬ 
terials when the ram pressure exceeds the gravi¬ 
tational restoring force per area. At the trunca¬ 
tion radius, the fraction of the ram press ure to the 
gravitational restoring force per area (ITakizawa 
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Fig. 6.- The thermal pressure profile of ”ID 1” centered on the mass center (open circles) and that of 
”ID 32” centered on the X-ray peak (open stars). The dotted (magenta) and dashed (green) lines indicate 
the ICM p ressure range of the sout hwest direction and azimuthal average excluding the southwest direction 
derived by Simionescu ct al. ( 2013l l. (A color version of this figure is available in the online journal.) 
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Fig. 7. - (a) The subhalo gas mass against weak-lening mass or hydrostatic mass. The open (red) circle 
and (blue) star marks indicate “ID 1” and “ID 32” , respectively. The error bars of the weak-lensing masses 
of the subhalos include the systematic errors. The filled (orange) circles and (cyan) stars are the gas mass 
vs. hydrostatic mass at the radius within which the mean density of the groups and clusters are the same 
overdensity as “ID 1” and “ID 32”, respectively. T he gas mass and h ydrostatic mass of the groups and 
clusters were calculated from Chandra’s results ( Vikhlinin ct al . m. (b) The same as left panel but 
temperature vs weak-lensing mass or hydrostatic mass. 
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u. r . 
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= 1.3 
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«e,ICM 


100 kpc/ V10 _o cm 

(1) 

is the useful parameter to investigate the present 
stripping effect. Here, the n e ,iCM, v ga i, r t , 
Msubhaio, and Afgas,subhalo are the electron density 
of the ICM, velocity of the subhalo, truncation 
radius, mass of subhalo, and gas mass of subhalo, 
respectively. By adopting those parameters of the 
subhalos and ICM, we calculated the ram pressure 
to the gravitational force per area ratio as a func¬ 
tion of the subhalo moving velocity and plotted in 
figure [8] We also estimated the ratio at the 0.6 r t , 
which corresponds to the border of the excess X- 
ray emission. The ram pressure should have been 
lower than the present value when infalling the 
outer regions, since the ICM density decreases 
with the distance from the cluster center. 

Since NGC 4807 and IC 4088 are located near 
the center of the mass contours of “ID 1” and 
“ID 32” (see figure [l]), respectively, these galaxies 
are likely representative galaxies of these subhalos. 
The recession velocities of NGC 4087 and IC 4088 
are 6989 km s" 1 and 7095 km s" 1 (NASA/IPAC 
Extragalactic Databas^l), respectively. Consider¬ 
ing that of the Coma cluster, 6925 km s _1 , the 
velocities in line of the sight of these galaxies are 
64 km s -1 and 170 km s _1 , respectively. As shown 
in figure [8j if these subhalos are moving toward 
the line of sight, the ratio of the ram pressure to 
the gravitational force in a unit area are orders of 
magnitudes lower than the unity. Even if these 
subhalos are moving with a inclination angle of 45 
degree, the ram pressure is not still effective to 
remove the gas of the subhalos. 

On the other hand, the infall velocities for 
”ID 1” and ”ID 32” are about 2000 kins” 1 , 
which is estimated using the best-fit NFW pro¬ 
file d erived from the w eak-lensing for the main 


-3 , 


halo ( Okabe et al. 2010ih This value is compa¬ 


rable with the infall velocity of the subcluster, 
NGC 48 39 group, whose infall v elocity is 1700±|gg 
km s -1 ( Colless fc Dunnlfl996 h. Adopting the in¬ 
fall velocity, at 1.0 r t , the ram pressure is higher 
than the gravitational force per area and enough 
to remove the gas, and at 0.6 r*, or at the bor- 


3 http://ned.ipac.caltech.edu/ 


der of the excess X-ray emission, the ram pressure 
is comparable to the gravitational force per area. 
AI^eaeisr ^3hal|6 ^4'fekmt^ ^y9 ,l ld ^oown the v/ky, 
lOhd/if^ljey a^e tQb^ijl^pe/penidfiSlJ&ikte Ahfc/line 
of sight with the infall velocity, their gas beyond 
0.6 r t can be stripped via ram pressure stripping. 

We also estimated the mass loss rate caused by 
Kelvin-Helmholtz instabilities since the ga s in the 
center of ”ID 32” seems to be removed. iNulsenl 


(1198211 estimated the mass loss rate caused by vis¬ 
cous stripping via Kelvin-Helmholtz instabilities, 
Mkh ~ Trr 2 D picMV. Here, and picm are the 
disk radius and gas density of the ICM, respec¬ 
tively. The mass loss rate can convert 


M, 


KH 


25 


( n « 'I 

( Td Vi 

f V \ 

V 10~ 5 cm~ 3 J 

V100 kpc J 

K 2500 km s _1 ) 


( 2 ) 

We assumed that rp is same as the truncation ra¬ 
dius. Using the velocity of line of sight of each 
subhalo, the results of the mass loss rate were the 
Mhk ~ 6 Mq yr -1 for “ID 1” and ~ 60 M 0 yr -1 
for “ID 32”. In contrast, adopting the infall ve¬ 
locity, the mass-loss rates of “ID 1” and “ID 32” 
are 180 and 680 M 0 yr _1 , respectively. Consider¬ 
ing the time scales for the mass-loss of 520 Myr 
and 390 Myr for ”ID 1” and ”ID 32”, respectively, 
infalling from the virial radius of the Coma clus¬ 
ter, the total mass loss from adopting the infall 
velocity the subhalo are about 9 x 1O 1O M 0 and 
3 x 10 n M e for “ID 1” and “ID 32”, respectively. 
These values of mass are more massive than cur¬ 
rent gas mass of the both subhalos. Therefore, de¬ 
stroying of gas by Kelvin-Helmholtz instability ex¬ 
plains the very low gas fraction of ”ID 1” and ”ID 
32” subhalos. However, the mass-loss rate would 
be smaller with magnetic fields which suppress the 
destruction by Kelvin-Helmholtz instability. 

4.3. The contribution of the X-ray emis¬ 
sion of subhalos to the ICM 

Suzaku enables us to measure the entropy 
profiles of the galaxy clus t ers out to the virial 
radius (George et al. 2009: Reiprich et al.1 2009t 


2010; Simionescu et al 


Bautz et al.ll2009: Hoshino et al.ll2010: Kawaharada et al 


20111; lUrban et al.1 l201lt 


Akamatsii£Wi]j2011 ,1 mnghre^et^I • 20.12: Walke r et al 


2012al; lAkamatsu et aLj 2012 Walkei^t al. 2012b; 

Sato et, all 20 12t I Walker etlil 20131: ISimionescu et, al 

2013: llchikawa et al. 2013t lUrban et ah 201 ll 
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Fig. 8.— The ratio of the ram pressure to gravitational force per area against the velocity of each subhalo. 
The solid and dashed lines are the ratio calculated at 1.0 n and 0.6 rt, respectively. The red and blue lines 
correspond to “ID 1” and “ID 32”, respectively. The vertical solid and dotted lines correspond to velocities 
of each subhalo assuming the 0° and 45° inclination angle, respectively. The infall velocities are plotted the 
red and blue dot-dashed lines for “ID 1” and “ID 32”, respectively. The horizontal solid line corresponds to 
unity. (A color version of this figure is available in the online journal.) 


Sato et al1l20l4[Okabe et al.l2014tJ : Mochizuki et al 


20151) . Contrary to the prediction of the ac¬ 
cretion shock heating model prediction, the en¬ 
tropy of galaxy clusters s how flat profiles beyond 
rsoo- ISimionescu et alJ (2011) interpreted that 
gas density in the outskirts are overestimated due 
to gas clumping and the entropy are underesti¬ 
mated. However, observing the Perseus cluster 
outskirts with Chandra , the number of detected 
sources is consisten t with the background sources 
( Urban et al •20 ll. With X-ray an d weak-lensing 
joint analysis, Okabe et al. ( 2014bl) discussed that 
entropy flattening of the outskirts of the galaxy 
clusters caused by the steepening of the tempera¬ 
ture profiles, rather than the flattening of the gas 
density. 


We studied the effect of subhalo luminosi¬ 
ties on the estimation of electron densities when 
these subhalos are not excluded from spectral 
analysis. The observed X-ray flux of the ex¬ 
cess emissions of ”ID 1” and ”ID 32” are (1- 
2)x 10 -13 erg s -1 cm -2 in the 0.5-2.0 keV energy 
band. Considering that the threshold of detection 
of point sources with a 10 ks Suzaku exposure 
is about 10 -13 erg s -1 cm -2 in the 2.0-10.0 keV 
range, if similar subhalos are located in other clus¬ 
ters, most of them would be below the detection 


threshold flux of Suzaku. We note that regions 
around subhalos like ”ID 9”, or the southwest 
subcluster around NGC 4839, can be easily ex¬ 
cluded from spectral analysis of clusters observed 
with Suzaku , because of their very high X-ray 
luminosities. 


Within projected distance 1.2-1.6 rsoo where 
the two subhalos are located, the ICM gas mass 
would be several times 10 13 Mq, which is esti¬ 
mated using the weighted average of the radial 
profi les of electron d e nsity observed with Suzaku 
by Simionescu ct al.l (201 ll) excluding the south¬ 
west direction, where the X-ray luminous subhalo, 
”ID 9” is located. Thus, the ICM gas mass is two 
orders of magnitude higher than the gas mass of 
excess emission of the "ID 1” and ”ID 32”. The 
0.5-2.0 keV luminosities of these two subhalos are 
a few times 10 41 erg s _1 and are negligible when 
comparing with the X-ray luminosity of the Coma 
cluster within 1.2-1.6 rsoo, ~ 10 43 erg s _1 , which 
is also estimated using the average of ICM tem- 
perature and normaliza tion observed with Suzaku 
( Simionescu et al ■'20 ll_ excluding the s outhwest 
direction. Although Okabe et al.l ( 2014a ) detected 
32 subhalos with weak-lensing observations, most 
of them are less massive and located within rsoo- 
As a result, we can conclude that X-ray emission 
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of subhalos does not affect the overestimating the 
gas mass of the Coma cluster. 

In order to evaluate the bias in the ICM temper¬ 
ature measurements, we extracted spectra again 
from FOVs observed around each subhalo, and fit¬ 
ted the spectra with ICM and background model. 
The temperature of the ICM did not change within 
statistic error range, since the total flux of subha¬ 
los is 1-2 orders of magnitude lower than that of 
the ICM. If number of subhalos are much higher, 
they may affect on the temperature and den¬ 
sity measurements. Therefore, we created the 
mock spectra assuming two thermal components 
for clumps and surrounding ICM emission, by 
changing the flux ratio of two components. When 
we simulated a sum of mock spectra of 2 and 1 keV 
emissions for the ICM and subhalo whose flux is 
half of that of ICM, respectively, the derived tem¬ 
perature with a single temperature model became 
1.2 keV, which is 40% lower than the original ICM 
temperature, and the electron density was also 
overestimated by 30%. This indicates the entropy 
to be low biased by 50%. Thus, if there are many 
clumps enough to increase the normalization of the 
ICM for changing the entropy profiles, the temper¬ 
ature bias is more significant. 

Simionescu et al.l ( 2013h derived the entropy 
of the Coma cluster beyond r R nn 4 7' ± 1') 


derived from Planck Collaboration et al. ( 2013 ). 


and the entropy profile is cons istent with the ac - 
cretion shock heating model ( Pratt et alJ l2010h . 


from IPlanck Collaboration ct alJ (12013 1. Since 
there is no evidence for the entropy flattening like 
other clusters and pressure excess, they suggested 
the gas clump are easily destroyed in a dynami¬ 
cal active cluster. Therefore, their discussion is 
consistent with our study that gas clump does not 
affect to the gas density estimation of the Coma 
cluster. 

5. SUMMARY AND CONCLUSIONS 

We observed the three massive subhalos, “ID 1”, 
“ID 2”, and “ID 32”, whi ch are detected from Sub¬ 
aru weak-lening analysis ( Okabe et al.ll2014a ) with 
Suzaku. The weak-lensing survey of subhalos in 
the outskirts of the galaxy cluster enable us to effi¬ 
ciently carry out follow-up X-ray observations gas 


subhalo candidates associated with weak-lensing 
detected subhalos. 

While the excess emission is seen around the 
center of “ID 1” mass contour, the “ID 32” sub¬ 
halo shows that the emission peak is shifted in the 
northern part or outer side from the Coma clus¬ 
ter center. In contrast to above two subhalos, the 
”ID 2” subhalo does not show any excess emission. 
The spectral analysis indicated that the tempera¬ 
ture of the subhalo gas is significantly lower than 
the surrounding ICM. By deprojecting the sur¬ 
face brightness profiles, we derived the gas mass 
of each subhalo. The total gas mass of “ID 1” and 
“ID 32” are 2 x 10 10 M 0 and 1 x 10 11 M 0 , re¬ 
spectively. Comparing with the kT — M to tai and 
M gas - Mtotai relation of regular galaxy groups, 
gas fractions of the subhalos are much lower than 
regular galaxy groups. Adopting the infall veloc¬ 
ity estimated from the best-fit NFW profile de- 
rived by weak-lensi ng analysis of the Coma cluster 
( Okabe et al.ll2010h ■ jeyond 0.6 times the trunca¬ 
tion radius of the subhalos, or at the border of 
the excess X-ray emission, the ram pressure is ef¬ 
fective to remove the gas. With the infall veloc¬ 
ity, total amount of the removed gas mass from 
the subhalos via Kelvin-Helmholtz instabilities are 
about 9 x 1O 1O M 0 and 3 x 1O U M 0 for “ID 1” and 
“ID 32”, respectively. The luminosities of subha¬ 
los are about two orders of magnitude lower than 
that of the Coma cluster outskirts and do not af¬ 
fect on the gas mass estimate of the ICM. 


Here, the r^nn derived from weak-lening analysis 


(Okabe et al. 2010) is well consistent with that 
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A. The CXB estimation 


For estimations of the Cosmic X-ray Background (C XB) level we extr acted the spectra beyond 110 


offset observations, which is corresponding to 2.5 r $oo ( Okabe et al. 2010h . excluding south region from 


the Coma cluster, whose observation logs are shown in table [2] We searched for point-like sources with 
“wavdetect” tool in CIAC0 in 0.5-2.0 and 2.0-5.0 keV images. We also excluded the area around the hot 
pixels @. The flux level of the faintest source was about lxlO -13 erg s -1 cm -2 in 2.0-10.0 keV with a 
power-law model of a fixed photon index, T =1.7. We assumed that the background emission were composed 
to two thermal Galactic emissions, the Local Hot Bubble (LHB) and the Milky Way Halo (MWH), and the 
Cosmic X-ra y background . The LHB and MWH were modeled non-absorbed and absorbed thermal plasma 
model (apec; ISrnith et~al 2001 ). The CXB was modeled an absorbed power-law model. We also convolved 
the Galactic ab sorption with ph otoe lectric absorption model, phabs. The column density was fixed to be 
8.5 x 10 19 cm 2 ( Ivalberla et al.lf2005h . Therefore, we modeled the spectra by the follow formula, constant x 
(apecLHB + phabs x (apecMWH + powerlaw)). Although the temperature of the LHB was fixed at 0.1 keV, 
the normalization was allowed to vary. The temperature and normalization of the MWH were free parameters. 
The photon index of the CXB was fixed at 1.4, and the normalization was allowed to v ary . The fitting result s 
of are summarized in table [5] The derived CXB normalization was cons istent with Kushino et al. (120021 ). 
The parameters of LHB and MWH are comparable with previous results (ISimionescu et al.ll2013l) . 


4 http: //cxc. harvard.edu/ciao/ 

5 http: //www. astro, isas.ac.jp/suzaku/doc/suzakumemo/suzakumemo-2010-01.pdf 


Table 4 

Suzaku OBSERVATION LOGS FOR THE CXB ESTIMATION. 


Field name 

Sequence No. 

Obs. date a 

(R.A., decl.) b 

J2000.0 

Exposure 0 

ksec 

East 110 / 

806037010 

2011-06-19T14:09:55 

13 h 07 m 48!6, 27°53'08"5 

11.1 

East 120 / 

806038010 

2011-06-20T01:23:51 

13 h 08 m 27?4, 27°53'06"0 

13.8 

NW 110' 

806045010 

2011-06-22T17:48:52 

12 h 56 m 18?6, 29°33'08"3 

14.1 

NW 120' 

806046010 

2011-06-23T06:03:44 

12 h 56 m 00?8, 29°4l'27('2 

11.9 


a Start date of observation, written in the DATE-OBS keyword of the event FITS files. 

b Average pointing direction of the XIS, written in the RA.NOM and DECJNTOM keywords of the 
event FITS files. 

c Exposure time after screening. 
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Table 5 

The background fitting results. 


Norm LHB 

fclMWH 

keV 

Norm^wH 

Norm CXB 

11 4 +1 ’ 4 

J " L 1.4 

q o-i +0.03 

0.841°;^ 

1 oi+ 002 
x ■ U_L —0.02 


a The normalization of the apec components divided 
the solid angle, Q u , assuming a uniform sky of 20 7 ra¬ 
dius, Norm = f n e niidV / [47t (1 + z) 2 D^} / Q u X 10 —17 
cm -5 4007r arcmin -2 . 

b The normalization of powerlaw is units of photons 
cm -2 s -1 keV -1 4007T arcmin -2 at 1 keV. 
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